Low-energy break in the spectrum of Galactic cosmic rays. 
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Measurements of the low energy spectrum of Galactic Cosmic Rays (GCR) by detectors at or 
near the Earth are affected by Solar modulation. To overcome this difficulty, we consider nearby 
molecular clouds as GCR detectors outside the Solar system. Using 7-ray observations of the clouds 
by the Fermi telescope we derive the spectrum of GCRs in the clouds from the observed 7-ray 
emission spectrum. We find that the GCR spectrum has a low energy break with the spectral slope 
hardening by Ar = 1.1 ±0.3 at an energy of E = 9±3 GeV. Detection of a low-energy break enables 
a measurement of GCR energy density in the interstellar space U = 0.9 ± 0.3 eV/cm 3 . 



Introduction. The spectrum of cosmic rays (CR) with 
energies E < 100 GeV measured by CR detectors inside 
the Solar System EH2 is not identical to the Galactic 
CR (GCR) spectrum due to the extinction of the low- 
energy CRs by the Solar wind [4j-[6| . Uncertainties in the 
knowledge of properties of the Solar wind, its termination 
shock, Heliosheath and Heliotail introduce uncertainties 
in our knowledge of GCR spectrum. 

The only possibility to measure the GCR spectrum 
unaffected by the Solar modulation is to consider nearby 
mass concentrations, like Giant Molecular Clouds (GMC) 
as natural CR detectors [TUTU] . The nearest GMCs form 
the Gould Belt, a ring-like structure of diameter ~ 1 kpc 
inclined at ~ 20° to the Galactic Plane QUE!]- CR in- 
duced 7-ray emission from the Gould Belt clouds was pre- 
viously detected by COS-B HUES], EGRET EH El] and 
by the Large Area Telescope (LAT) on board of Fermi 
satellite PIES]- 

The bulk of 7-ray emission from the GMCs is produced 
by GCR interactions. The GCR spectrum may be re- 
constructed using the measured 7-ray spectrum from the 
clouds combined with the known relevant particle physics 
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TABLE I: High Galactic latitude Gould Belt clouds consid- 
ered in the analysis .l s , b s and O s are the coordinates of the 
center and radius of regions from which source counts are 
collected, lb, bb are the coordinates of the centers of 5° re- 
gions used for background estimation. Masses M (in units of 
10 5 Mq) and distances D (in parsecs) to the clouds are from 
Ref. [IT]. 



of pion production and decay [2TJI |2"2] 1 . In what follows 
we report such a measurement. 

Data selection and analysis. Clouds from the Gould 
Belt span large angular sizes 6 ~ 1° — 10°. Most of 
these clouds lie in the Galactic Plane, so that the 7-ray 
emission is superimposed on the diffuse Galactic emis- 
sion. This makes the analysis of characteristics of cloud 
7-ray emission difficult: the diffuse Galactic emission is 
variable on the angular scales comparable to and has 
similar spectral characteristics. Clouds at high Galactic 
latitudes are separated from the Galactic Plane. Analy- 
sis for these clouds can thus be done in a straightforward 
way. Taking this into account, we concentrate on the 
study of high Galactic latitude (|6| > 10°) clouds listed 
in Table H] 

We use the LAT data collected between August 4. 
2008 and July 15, 2011. We filter the data with Fermi 
Science Tools 2 (software version v9r23pl and data se- 
lection p6_vll) using gtselect, gtmktime and gtbin. We 
retain only 7-ray events (dataclean events) at zenith 
angle < 100°. We use the aperture photometry method 
for the spectral analysis, by collecting the events from 
"source" regions and comparing the total number of 
counts in each source region with the number of back- 
ground events, estimated from nearby "background" re- 
gions at the same Galactic latitude. The lists of "source" 
and "background" regions is given in Table [I] The expo- 
sure is calculated using gtexposure tool. 

Large regions occupied by the clouds contain point 
sources. Emission from the point sources is superimposed 
on the diffuse emission from the clouds. To subtract the 
point sources, we use the list of sources from the two year 
exposure of LAT [23] . 

The point spread function (PSF) of LAT becomes com- 
parable to the size of the clouds at energies < 0.3 GeV for 
photons pair converted in the "front" layer of the LAT 
and at < 0.7 GeV for photons converted in the "back" 



1 A comparison of different codes for calculation of pion production 
in pp interactions was recently done in Ref. |21| 

2 http: / / fermi.gsfc.nasa.gov/ssc / data/ analysis /scitools / 
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FIG. 1: Spectrum of diffuse 7-ray emission from high Galactic 
latitude GMCs. Red thin curve is the spectrum of 7-ray emis- 
sion calculated assuming powerlaw CR spectrum. Orange 
hatched region is "local atomic hydrogen" component from 
the Ref. [23], renormalized by a factor 23.7 to match the nor- 
malization of the GMC average spectrum. Grey shaded area 
shows model spectra calculated assuming broken powerlaw 
GCR spectrum. 

layer. We use only front photons for the analysis between 
0.3 and 0.7 GeV. We do not use the E < 300 MeV data. 
Gamma-ray spectrum of the clouds. All the clouds 
listed in Table |T] are detected as extended sources with 
LAT. The morphology of 7-ray emission follows the mor- 
phology of CO emission as inferred from the CO maps 
[TT1 IT3] . The 7-ray flux is proportional to the CO 
integrated intensity as expected if the flux scales as 
F ~ M/D 2 with M and D being the cloud mass and 
distance [7HTD]. Images and spectra of individual clouds 
can be found in the Supplemental Material. 

We verified that the spectra of individual clouds are 
consistent with each other and with the " local atomic hy- 
drogen" component of diffuse Galactic emission derived 
in Ref. [23] (Fig. [I]). This implies that all the high- 
latitude Gould Belt clouds are "passive" CR detectors, 
with no on-going particle acceleration and no modifica- 
tion of CR propagation inside the clouds (i.e. the CR 
diffusion coefficient does not change significantly in the 
clouds [9]). The average spectrum of the clouds is shown 
in Fig. @ 

From Fig. [T] we see that the 7-ray spectrum exhibits 
a break at ~ 2 GeV. Several explanations of this break 
can be considered: a break in the spectrum of CR nu- 
clei or CR electrons/positrons or a break in the 7-ray 
production cross-section. 

7-rays with energies ~ 2 GeV are produced by pro- 




FIG. 2: 68% confidence ranges for -EWeak, Ti,r2 and a pa- 
rameters. Black dot shows the best-fit. Light-blue shaded 
region is for a — > 00, grey-shaded region is for finite a. 

tons with energies much higher than the pion produc- 
tion threshold, where proton-proton interaction cross- 
section grows logarithmically with energy. The break 
could also not be related to a feature in the electron 
Bremsstrahlung cross-section since this also grows log- 
arithmically at GeV energies. The most reasonable ex- 
planation for the ~ 2 GeV feature is that it is related to 
a feature in the CR spectrum. 

CR spectrum. We have reconstructed the spectrum 
of CRs from the 7-ray spectrum using the parametrized 
pion production spectra calculated in Ref. 20J . Fig. 
[T] shows a comparison of the observed 7-ray spectrum 
with that produced by a powerlaw distribution of CR 
protons/nuclei. The spectrum produced by a powerlaw 
CR spectrum has a peak at ~ 300 MeV and the model 
over-predicts the 7-ray emission below ~ 1 GeV. The 
sub-GeV flux could be suppressed if the GCR spectrum 
hardens below ~ 10 GeV. 

To find the details of the low-energy hardening we 
model the spectrum as a broken powerlaw dNcn/dE ~ 

(E/E Bl ) ri / (l + (E/E hle ^) a f 2+ri)/ ° with low and 

high energy slopes T%, T2, break energy inbreak and sharp- 
ness of the break a. Assuming negligible Bremsstrahlung 
contribution, a satisfactory fit to the 7-ray spectrum is 
found, with x 2 /d.o.f = 23/21 which corresponds to a 34% 
probability for model to be the proper description of the 
data. The range of model 7-ray spectra consistent (at 
68% confidence level) with the GMC data is shown in 
Fig. [11 

The conclusion about the presence of a low-energy 
break in the CR spectrum is not altered if a non- 
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negligible Bremsstrahlung component is considered. The 
maximal Bremsstrahlung component is produced when 
electrons (primary CR electrons and secondary electrons 
produced in interactions of CR nuclei) loose all their 
energy via Bremsstrahlung before leaving the clouds. 
This could happen if electron diffusion through the 
clouds is slower than in the ISM. This is not real- 
ized in the Gould Belt, but considering the maximal 
Bremsstrahlung component allows an estimate of influ- 
ence of Bremsstrahlung-related uncertainties. Includ- 
ing the maximal Bremsstrahlung component we find a 
satisfactory fit (x 2 /d.o.f ~ 24/21) which is achieved 
with a sharper break in the CR spectrum, because 
Bremsstrahlung contributes mostly to the lower energy 
part of the spectrum at E < 300 MeV [S]. 

To obtain the estimates for inbreak, Ti, L2, cr we follow 
Ref. [25.. Projections of the 68% confidence intervals 
for the model parameters onto (£Weak, c), (-EWcak, I^), 
(ri,r 2 ) and inbreak, cr planes are shown in Fig. [2] The 
sharpness of the break, er, is constrained from below. 
The best-fit model is the model of the form dNcn/dE ~ 
(E- T \ E < E Br ) U {E- t \ E > E Br ). The confidence 
regions for model parameters in the " sharp break" model 
are shown as light-blue shaded regions in Fig. [2j Project- 
ing these regions onto coordinate axis one finds -Ebrcak = 
9±3 GeV and T 2 = 3.03 ±0.17, I\ = 1.9±0.2. Consider- 
ing the model with maximal Bremsstrahlung component 
one finds a reduced estimate for Ti = 1.7±0.2. Leav- 
ing parameter a free increases the errorbars to -Ebrcak = 



9±1 GeV, T 2 = 3.03 
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The PAMELA collaboration has reported a measure- 
ment of the CR spectrum between 50 GeV < E < 
200 GeV with the spectral slope of T ~ 2.85 [2]. LAT 
measurements are mostly sensitive to the part of the CR 
spectrum below 200 GeV. The slope of the CR spectrum 
L 2 is consistent with the PAMELA measurement in the 
50 - 200 GeV range. Below 50 GeV the CR spectrum 
derived from Fermi measurement is still consistent with 
a powerlaw down to ~ 10 GeV, while the spectrum mea- 
sured by the PAMELA deviates from the powerlaw. 

Contrary to the shape, the normalization of the GCR 
spectrum is more difficult to measure using the 7-ray 
data. The problem here being the large uncertainty (by 
a factor of ~ 2) of the of amount of target material in 
the clouds. To remove this uncertainty, we normalize the 
CR flux above 200 GeV in the PAMELA data. This is 
justified because the effect of the Heliospheric distortion 
at this energy is minor. The GCR spectrum normalized 
in this way is shown by dark shaded curve in Fig. [3] 

An alternative possibility to find the normalization 
of GCR spectrum is to infer the 7-ray emissivity per 
hydrogen atom which is determined by the CR den- 
sity. Such a measurement relies on an estimate of 
the hydrogen column density from the CO maps us- 
ing the CO-to-N# 2 conversion "X- factor". Assuming 
X = 1.8 x 10 20 cm _2 K _1 krn _1 s 13] we derive the nor- 
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FIG. 3: Light grey: GCR spectrum inferred from the 
LAT observations. Dark grey: GCR spectrum inferred from 
the LAT observations with normalization fixed at 200 GeV 
to PAMELA. Thin data points are proton (light blue), he- 
lium (cyan) and the total (blue) CR spectrum measured by 
PAMELA 0. Olive data points in 0.1-0.35 GeV range are 
Voyager data for CR flux beyond the Solar wind termination 
shock [3]. Solid curve is GCR spectrum reconstructed from 
Voyager [31 [BJ. 

malization of the CR spectrum shown by the light grey 
band in Fig. [3] A constant nuclear enhancement factor, 
re ~ 1.5 — 1.8, which accounts for the nuclear composi- 
tion of the CR flux and of the interstellar medium, was 
assumed [551 EH]- Statistical uncertainty of such mea- 
surement (~ 10%) is much smaller than the systematic 
uncertainty related to uncertainty of the X-factor [19] 
(~ 40%) and by the uncertainty of re (~ 20%) J2HE01)- 
Discussion. The LAT observation of the high Galactic 
latitude clouds from the Gould Belt shows that the steep- 
ening of GCR spectrum below ~ 200 GeV persists down 
to - 10 GeV. The PAMELA spectrum below - 50 GeV 
deviates from the GCR spectrum derived from the LAT 
data. This could be attributed to the distortion of the 
GCR spectrum in the Heliosphere. In the conventional 
modeling, the Heliosphere is assumed to affect the CR 
flux only below ~ 10 GeV, in the energy band where the 
Solar modulation, or time variability of the flux is ob- 
served pQ. However, the heliospheric effects might affect 
the CR spectrum up to TeV energies at which the gy- 
roradius of CRs becomes comparable to the size of the 
Heliosphere, ~ 100 AU. At 50 GeV, the gyroradius is 
R L = E/eB ~ 1 [E/50 GeV] [5/10 ^Gf 1 AU, compara- 
ble to the size of magnetic structures in the outer Helio- 
sphere recently revealed by the Voyager spacecrafts [31] . 
The possibility that the Heliosphere influences the CR 



4 



flux at the energies higher than 10 GeV, which though 
physically justifiable, needs to be further investigated. 
Residual influence of the Heliosphere might be present 
up to still higher energies, ~ 1 — 10 TeV, at which the 
anisotropy of the CR flux in the direction toward and 
opposite of the Heliotail is observed [55] . 

The GCR spectrum breaks by Ar ~ 1 below ~ 
10 GeV. Such a break was not reported before, although 
indications for the existence of a break, based on compar- 
ison of the CR data with GCR propagation models, were 
discussed [STJ[2E1[1I] ■ ^ possibility of existence of a break 
in the spectrum of CR electrons in the same energy range 
was recently discussed [23] . The main difference between 
our result and that of Refs. [21] [26] [27] is that we ob- 
tain a direct and model-independent measurement of the 
break, without any assumptions on the distribution of 
sources in the Galaxy and details of CR propagation. A 
potential uncertainty of our measurement is possible ef- 
fect of CR propagation in the clouds on the measured 
spectrum. Absence of cloud-to-cloud variations of the 
position of the break argues against such a possibility. 
Small cloud-to cloud variations might be detectable with 
deeper exposure by LAT. 

The hard slope of GCR spectrum below the break, 
Ti < 2, is important because it insures a finite en- 
ergy density of GCRs [5|. Upper limits on the GCR 
flux in the 0.1-0.35 GeV range derived from the Voy- 
ager measurements [3J [S] (Fig. [3]). Combining the 
Voyager constraint with the LAT measurements we de- 
rive a measurement of the GCR energy density, Ucr — 
0.9 ± 0.3 eV/cm 3 . Such an energy density is in equipar- 
tition with the energy density of magnetic fields Ub — 
1 [B/6 /iG] 1 / 2 eV/cm 3 and turbulent motions U tur b — 
1 [n/SAf/1 cm" 3 ] [vturb/20 km/sf eV/cm 3 of the ISM 
with density ujsm an d turbulent velocity v tur b, a fact 
which points to a physical coupling between the three 
ISM components. 

The detection of a low-energy break introduces a new 
energy scale into the CR physics. Taking into account 
that the new energy scale is not far from the "natural" 
scale of proton rest energy m p c 2 , we have tested a pos- 
sibility that the new scale could be reduced to m p c 2 if 
the broken powerlaw model of CR spectrum is changed 
to an alternative model of the form dNcR,/dE ~ f3 Pl lZ P2 , 
where /3 and 1Z are CR velocity and rigidity. Such shape 
of GCR spectrum does not have additional scale except 
for TOpC 2 . We find that such a GCR spectrum is inconsis- 
tent at > 4a with the LAT data, with the best-fit value 
X 2 /d.o.f. > 3 for 23 d.o.f. 

A change of the CR spectral slope might be related 
to the physics of CR sources (e.g. characteristic maxi- 
mal/minimal/break energy of CRs produced by a source 
population J3T] HS1 HZ] , or injected by annihilation/decay 
of Dark Matter) or to a change of the propagation for 
CRs (e.g. transition from convective to diffusive regime 



[H [5S] or change in the energy dependence of the dif- 
fusion coefficient [26, 32J and/or diffusive re-acceleration 
possibly combined with an intrinsic break in the source 
spectra [2"Tl|g2]b 

A low-energy cut-off in the source spectra could occur 
e.g. if the lower-energy particles are trapped by mag- 
netic fields inside the sources. CRs of energies ~ 10 GeV 
propagate through the Galaxy in a diffusive way by scat- 
tering on turbulent inhomogeneities in ISM. A feature in 
the turbulence spectrum at a length scale At might pro- 
duce a break in the CR spectrum at an energy at which 
the Larmor radius Rl = Ecn/eB (B is magnetic field 
in the ISM) is comparable to Ay. Scattering and/or ab- 
sorption of the lower energy CRs would be determined by 
the energy-independent geometrical cross-section of the 
smallest ISM inhomogeneities. In such a model, measure- 
ment of the break energy -Bbroak ~ 10 GeV implies the de- 
tection of a feature in the distribution of inhomogeneities 
of ISM at the length scale X T ~ E break /eB ~ 1 AU. 

Suppression of the low energy CR flux might also oc- 
cur through efficient CRs interactions with the ISM on 
the time scale of proton-proton interactions t pp ~ 3 x 
10 7 [njg M /l enr 3 ] 1 yr. During this time, CRs could 
spread over a region of the size R(E) < yjD(E)t pp ~ 

1 [E/l GeV] ' 3 [n ISM /l cm" 3 ] -172 kpc around a CR 
source (assuming D(E = 1 GeV) ~ 10 28 cm 2 /s for the 
CR diffusion coefficient). If the distance to the nearest 
GCR accelerator is in the kiloparsec range, energy losses 
would efficiently remove CRs with energies below several 
GeV from the locally observable GCR flux. This mech- 
anism of suppression of the GCR flux could work only 
if the last episode of injection of GCR within a 1 kpc 
volume occurred not later than tccR ~ tpp ~ 3 x 10 7 yr 
ago. Remarkably, this estimate is close to the age of the 
Gould Belt, which was formed in an explosive event some 
t GB ~ 3 x 10 7 yr ago [H]. 
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Supplemental Material 



SI. Imaging analysis 



All the clouds listed in Table 1 are detected as extended sources with LAT. Comparison of the countmaps shown in 
Fig. S1-S7 for the energy band E > 1 GeV with the CO intensity [11,13] shown by contours reveals a good correlation 
between the 7-ray and CO emission. Regions used to estimate the source and background fluxes for each cloud (Table 
1) are shown as white solid and dashed circles in Figs. S1-S7. 



I 2 3 4 5 B 7 _ S 1 

Figure SI. Count map of Perseus OB2 cloud region in the energy range E > 1 GeV smoothed with a Gaussian of the 
width 0.3°. Green contours show CO emission intensity with the levels 5, 15, 25, 35 and 45 K km/s. Crosses show 
positions of sources from the two-year Fermi catalog [24] . Solid and dashed white circles are regions used to estimate 
the source and background fluxes. 



G 




Z 3 4 5B7H3 



Figure S2. Same as in Fig. SI but for the Taurus cloud. 
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Figure S3. Same as in Fig. SI but for the Orion cloud region. 
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Figure S4. Same as in Fig. SI but for the Chameleon cloud region. 
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Figure S5. Same as in Fig. SI but for the p Ophiuchus cloud region. 



Figure S6. Same as in Fig. SI but for the r Cr A cloud region. 
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Figure S7. Same as in Fig. SI but for the Cepheus cloud region. 
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S2. Spectral analysis 

S2.1 Perseus 0B2 



i 1 — i 1 — i 
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Figure S8. Spectra of Perseus 0B2 (left) and Taurus (right) clouds. Black data points show the total background 
subtracted spectra of the source regions, red data points are the spectra of diffuse emission. Other data points 
are spectra of point sources from the two-year Fermi catalog [24] which are subtracted from the total source region 
spectrum in the estimate of diffuse emission spectrum. Left panel: blue and green thin data points show the spectra 
of NRAO 140 and 1RXS J033348. 7+29164. Right: green: BZB J0433+2905; blue: CRATES J0456+2702; magenta: 
BZB J0440+2750. Points without y errorbars are upper limits. 

From Fig. SI we see that two point sources, identified with a quasar NRAO 140 and an X-ray source 1RXS 
J033348. 7+29164 [24], are situated in the middle of diffuse emission region. Spectra of these sources reported in the 
two-year Fermi catalog are shown by the blue and green thin data points in Fig. S8, left. We subtract fluxes of these 
sources from the overall flux of the source region to obtain the estimate of diffuse emission from Perseus OB2 cloud 
(shown in red in Fig. S8, left panel). 



S2.2 Taurus cloud 

Taurus molecular cloud is closest to the Solar system GMC situated in the center of the Gould Belt [11]. It spans 
a large angular size of about 15° x 15° in the direction close to the Galactic Anticenter. A large sub-cloud situated 
in the Auriga sky region is at Galactic latitude \b\ < 10°. We do not consider this part of the cloud in our analysis. 

An unidentified source 2FGL J0440.5+2554c is clearly associated to the cloud, being positionally coincident with 
the excess of CO emission (Fig. S2). We consider this source as a part of diffuse emission from the cloud. 

The source region of Taurus cloud, used for spectral extraction, contains two blazars BZB J0433+2905, BZB 
J0440+2750 with the spectra shown in Fig. S8. Spectra of these blazars are subtracted from the overall spectrum of 
the region in our spectral extraction procedure. 



S2.3 Orion A.B and Mon R2 clouds 

The Orion molecular cloud is conventionally divided onto A and B sub-clouds, which are clearly distinguishable in 
the count map shown in Fig. S3. 



10 



E 




E [HeVj 



Figure S9. Same as in Fig. S8, but for Orion A (left) and Orion B (right) clouds. Left panel: blue and green are 
spectra of Crates J0538-05231; green: BZQ J0529-0519. 

The unidentified Fermi sources 2FGL J0547.1+0020c, 2FGL J0547.5-0141c, 2FGL J0543.2-0120c, 2FGL J0541.8- 
0203c, 2FGL J0534.9-0450c, 2FGL J0534.8-0548c are most probably associated with emission from the Orion clouds 
themselves, as can be judged from the correlation of source positions with the excesses in the CO emission (Fig. S3). 
We do not subtract the spectra of these sources from the overall emission spectrum from the Orion cloud source 
regions when estimating the spectrum of diffuse emission. 

Sources BZQ J0529-0519 and CRATES J053812-05231 situated inside the region used to estimate the source flux 
are identified with Flat Spectrum Radio Quasars (FSRQ) in the two-year Fermi catalog [24] are subtracted from the 
overall source region spectrum (Fig. S9). 

Mon R2 cloud region contains a point source OH-10, which is clearly not associated to the cloud emission and is 
subtracted in the spectral extraction procedure (Fig. S10). The cloud has small angular extent (2° radius circle is 
used in the spectral extraction). Taking this into account we consider only the data at the energies where the LAT 
PSF becomes smaller than the size of the region used for the spectral extraction and disregard the data below 1 GeV 
for this cloud (Fig. S10). 



S2.4 Chameleon cloud 
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Figure S10. Same as in Fig. S8 but for Mon R2 (left) and Chameleon (right) clouds. Left panel: green data are the 
spectrum of OH -10. Right: green data are the spectrum of PKS 1057-79, blue data are for PKS 1221-82. 
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Contrary to Perseus OB2, Orion A, B and Mon R2 clouds, which are situated in a low background Galactic 
Anticenter region, the Chameleon cloud is situated much closer to the Galactic Center, at Galactic longitude I ~ 300°. 
However, at Galactic latitude b ~ —15° the diffuse 7-ray emission from the cloud region is clearly separated from the 
Galactic Plane diffuse emission as it is clear from Fig. S4. 

Background sources PKS 1057-79, PKS 1221-82 and PKS 1133-739 are clearly not associated to the cloud. We 
subtract the spectra of these sources from the overall source region emission when estimating the spectrum of diffuse 
emission from the cloud (Fig. S10). 



S2.5 Rho Ophiuchus cloud 

p Ophiuchus cloud is situated in the Galactic Center region with strongly variable Galactic diffuse emission extending 
to high Galactic latitudes and numerous point sources (Fig. S5). 

A pulsar PSR J1614-2230 could be clearly identified as a point source unrelated to the cloud diffuse emission. PKS 
1622-253, a source identified with a blazar, not associated to an excess of CO emission, is, most probably, also a 
background source. However, the spectrum of this source, shown in Fig. Sll, is similar to the spectrum of diffuse 
emission from the Gould Belt clouds. This might indicate that the source spectrum reported in the two-year Fermi 
catalog [24] is contaminated by the diffuse emission from p Ophiuchus cloud. In this case subtracting the source 
spectrum from the overall spectrum extracted from the p Ophiuchus region we slightly under-estimate the diffuse flux 
from the p Ophiuchus cloud by up to 20%. 

Situation is more complicated with CRATES J1626-7638. This source is identified in the two-year Fermi catalog 
as a blazar. At the same time, it is spatially coincident with a strong excesses in the CO map. Moreover, the source 
spectrum repeats the spectrum of diffuse emission from the Gould Belt clouds (Fig. Sll). Taking this into account 
this uncertainty, we do not subtract the spectrum of CRATES J1626-7638 from the overall spectrum of p Ophiuchus 
source region when extracting the spectrum of diffuse emission from the cloud. 
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Figure Sll. Left: spectrum of p Ophiuchus cloud (red). Green: spectrum of PKS 1622-253; magenta: spectrum of 
PSR J1614-2230; blue: spectrum of CRATES J1626-7638. Right: spectrum of r Cr A cloud. 



S2.6 Corona Australis cloud 

Corona Australis (r Cr A) cloud is located in the Galactic Center region at Galactic latitude b ~ —20°. In spite of 
the low mass of the cloud, 7-ray emission is detectable because of the proximity of the source. 

Maximum of the 7-ray signal coincides with an excess in CO emission (Fig. S6) and is located at the position of 
an unidentified Fermi source 2FGL J1904.9-3720c. We propose that this source is associated with the cloud emission 
and do not subtract the spectrum of this source from the overall spectrum of the source from the total source region 
spectrum (Fig. Sll, right). 
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S2. 7 Cepheus cloud 
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Figure S12. Spectrum of Cepheus cloud (red). Grey: average cloud spectrum. 

Cepheus cloud is a bright and very extended 7-ray source situated just above the Galactic Plane in vicinity of 
Cygnus regon (Fig. S7). An unidentified source from the two-year Fermi catalog [24] is clearly associated to a strong 
excess in the CO emission map and we propose that this source should be identified as a clump in the cloud. We do 
not subtract the spectrum of this point source from the in our estimate of diffuse emission from Cepheus cloud (Fig. 
S12). 



